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Abstract. The vortex-lattice melting transitions in two typical iron-based high-Tc 
superconductor Ba{Fe\- x Co x )2As2 (122-type) &ndNd(Oi- x F x )FeAs (1111-type) for 
magnetic fields both parallel and perpendicular to the anisotropy axis are studied 
within the elastic theory. Using the parameters from experiments, the vortex-lattice 
melting lines in the H-T diagram are located systematically by various groups of 
Lindemann numbers. It is observed that the theoretical results for the vortex melting 
on both Ba(Fe\- x Co x )iAs2 and Nd(0\^ x F x )FeAs for parallel fields agree well 
with the recent experimental data. The future experimental results for the vortex 
melting can be compared with the present theoretical prediction by tuning reasonable 
Lindemann numbers. 
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1. Introduction 

Recently, the newly discovered iron-based superconductors have attracted considerable 
scientific interest both experimentally Pfl |3l El SI ISl El IZl IHl El 1101 IIH HSl H31 H31 HSl HSl 11^ 
and theoretically [HI HH EDI EH E21 EH]. First, the novel superconductivity in 1111 
phase FeAs superconductors was reported experimentally, giving a new path to high 
temperature superconductivity. The LaOFeAs under doping with F~ irons at the 
2 ~ sites was found to exhibit superconductivity with T c = 26K [TJ, and the 
superconductivity was also observed with holes doping [2]. Then T c was surprisingly 
increased up to above 40K when La in LaO\- x F x FeAs was substituted by other 
rare earth elements [31 H]. It was quickly observed later that T c is about 55K in 
SmOi_ x F x FeAs [51 [6] and Gd\_ x Th x OFeAs [7]. It is a first non-copper-based 
superconductors in which the maximum critical temperature is much higher than 
the theoretical value predicted from BCS theory [21]. On the other hand, the 122 
phase iron-based superconductor BaFe2As2 was discovered more recently [8], and 
the superconducting critical temperature was found to be as high as 38K by hole- 
doping. It is observed that the electron doping of which by Co [9] and Ni [10] 
also induce superconductivity. Although there are some theoretical studies on iron- 
based superconductors, the type-II superconductivity as well as the mechanism of 
supeconductivity are not well understood to date. 

The vortex-lattice solid (glass with random pinning) state with zero linear resistivity 
is crucial for the application of high-T c superconductors, thus the melting of vortex- 
lattice in bulk type-II superconductors is of great significance [23, [213 123 I2S I2S]. The 
main aspect of Lindemann criterion suggests that the lattice melts when the root mean 
square thermal displacements of the components of a lattice reach a certain fraction of 
the equilibrium lattice spacing, such a criterion was first adopted to study the vortex- 
lattice melting transition in type-II superconductors with a magnetic field parallel to 
the anisotropic axis, then this approach was used to draw the melting lines in the case 
of magnetic field perpendicular to the anisotropic axis [30j [3TJ [32l [33] . Since the upper 
critical field is also high [131 H2] in iron-based high-Tc superconductors, the thermal 
fluctuation may drive the vortex-lattice to a vortex liquid in a field far below the upper 
critical field [30] through vortex melting. We will extend the elastic theory to study the 
vortex melting in these newly discovered superconductors. 

In this paper, using the parameters measured in recent experiments, we study 
the vortex-lattice melting transitions for two typical iron-based layered superconductors 
Ba(Fe 1 - x Co x ) 2 As 2 (122-type) with low anisotropy PH2] and Nd{0 X - x F x )FeAs (1111- 
type) with high anisotropy [HEE], in the framework of the elastic theory. The melting 
lines for magnetic fields both parallel and perpendicular to the anisotropic axis are 
systematically located with different groups of Lindemann numbers. A comparison 
with current experimental findings is made. The present paper is organized as follows. 
In Section 2, we introduce the theoretical method used in this work, Section 3 presents 
the main results, finally, we give a short summary. 
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2. Elastic theory 

We will consider the field B both parallel and perpendicular to the anisotropy axis (i.e. 
c-axis in this paper), the elastic theories in both cases are presented respectively. 

2. 1 . Thermal fluctuations 

Whether the field is parallel or perpendicular to c-axis, for ideal triangular vortex-line 
lattice, the free energy in elastic theory can be expressed in an unified way with quadratic 
terms of the deviation vector u = (u x ,u y ) describing the fluctuations of vortices from 
their equilibrium positions [251 1251 EH EH 1221 132] 

l r d 3 k 

-u-Cu, (1) 



2 J (2tt) 3 ' 

The matrix C for fields parallel to the c-axis is different from that for fields perpendicular 
to it. We denote C c and C ab to be the elastic matrix for the fields parallel and 
perpendicular to c-axis, respectively, which are given as follows 

C c _ ( c lK + c m k\ + c 44 k 2 z c L k x k y \ 

V c L k x k y c L k 2 y + c m k\ + cukl J ~ ' 



and 

/ , b 2 _L r h . b 2 -L r h b 2 r-,,bb \ 

(3) 



iab _ ( c uk 2 + Cggfcy + c\ 4 k 2 z c\\k x k y 

\ C\\k x ky c mk 2 x + c u ky + c e i4 k 2 z 



In matrix C c , k\ = k 2 + k 2 , cqq,cl, and C44 are the wave-vector-dependent shear, 
bulk, and tilt elastic moduli [26] [Ml ES], respectively, which are determined as follows 

B 2 M 1-b 1 
c 44 (k) - ^ijfj 26^-ljfe2 + (M/M z )(k 2 + m 2 ) + 1J ( } 



and 



5 2 l-6 r k 2 + {M/M z )m\ 



4tt IbK 2 1 (k 2 + m\) [k\ + {M/M z ) (k 2 + m 2 )} 

~ k 2 ± + {M/M z )k 2 + m 2 ^ (5) 

Where we have defined m\ = and m| = M is a quasi-particle effective mass in 
x?/-plane, M z describes along c-axis, b = B/B c2 with B c2 is the upper critical field, and 
k = X±/k±. The bulk modulus is cl = en — cqq with shear modulus 

B 2 c2 b(l-b) 2 

The matrix C ab also contains some elastic moduli. Since the anisotropy exists here, 
tilt and shear modulus are not isotropic any more. For example, c^k) {c%J)s)) is tilt 
modulus along (perpendicular to) c-axis. Similarly, Cg 6 (c| 6 ) represents shear modulus 
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parallel to (perpendicular to) the anisotropy axis c. For the detailed expressions for 
elastic moduli, one may refer to Ref. [33] 

The thermal fluctuations of the vortices are given by inverting the kernel C(k) as 
follows 

«) = / dk[Ctr\k), a = x,y,(3 = c,ab (7) 

The integrations in Eq. (7) are over < k z < oo and within the first Brillouin zone 
(BZ) for k x and k y . To be specific, we consider a lattice structure as the low temperature 
equilibrium state as shown in Fig. 1 of Ref. [33], where the lattice spacings are a x = a 



and a y = aJ3/2 for fields parallel to c-axis, a x = aj ^ with x \\ c and a y = 



for fields perpendicular to c-axis with 7 2 = m c /m ab and a = y 2$o/v / 3-£>. 

For parallel fields, we consider the mean-square displacement of a vortex lattice 
from the equilibrium d 2 (T) = u 2 + u 2 which can be written as 

rf2(T) = jinf I dk[ c m ki + c u ki + Cll ki + c u ki ] (8) 

For convenience, we introduce the dimensionless wave vector q = (k x /A x , k y /A y , k z /A), 
where A x>y are the wave numbers at the edges of the first BZ. Explicitly, they are given 
as A x = Air/3a and A y = 27r/a\/3 for fields parallel to c-axis and A x = 47r v /7~/3a and 
= 27r/ay / 37 for fields perpendicular to c-axis. The unit of wave number in k z di- 



rection is taken as A = y4irB/§ = ^87r/y/3a 2 . The thermal fluctuations are then 
expressed as 

W-^JttZ^-Si^I* forBWc (9) 

and 

^ = jpT^s J ^^(q) = ^j0^ Iaa ' a = x ^y^ f° rB ±c ( 10 ) 

with the normalized matrix C(q) = C(q)/-B 2 A x A r 
2.2. Lindemann criterion 

The Lindemann criterion presumes that the lattice melts, when the root mean square 
thermal displacements of the components of a lattice reach some fraction of the 
equilibrium lattice spacing. For two cases, we write the usual isotropic Lindemann 
criterion for parallel fields, 

(d 2 ) = c 2 a 2 . (11) 

and the anisotropic one for perpendicular fields, 

(u 2 x ) = c 2 x a 2 x , (ul) = c 2 y a 2 y . (12) 
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with c x and c y are two Lindemann numbers for two transverse directions. Combining 
the Lindemann criterion and the elastic theory, we can get the melting equations for 
two cases, 



7- 

and 



2 "~" " ^ " forB\\c (13) 



a/1 — t \/3£d Idd 



= ( — ) c~ — = — - — , a = x,y, jorB _L c (14) 
a/1 — t a a V3e V 

with e the Ginzburg parameter e d = 167r 3 /t 4 (/c B T c ) 2 / ( l ) QiJ^ 2 (0) and e = 
16n 3 K\k B T c ) 2 /<S> 3 H$(0) . 

2. 3. Layer pinning effect 

When fields perpendicular to c-axis, the effect of layer pinning reduces fluctuations in 
both directions and induces an additional momentum-independent term to the elastic 
matrix in Eq. (3) such that 

C ab = f Cl1 ^' C 66^y + C 44^z + © C\ik x k y \ 

p \ C\\k x ky C\\ky + Cggfc 2 + C44A; 2 y 

where 

-®- = V^Aaj^ 1 e -A(4) a (16) 
A 2 fi 2 /3 A « 2 7 s & 2 (l-t) 3/2 1 J 

is proportional to the critical depinning current [32], with = ^ C (T = 0), ~ 1.16 
and s the layer separation. 



3. Results 

3.1. Vortex melting in Ba{Fe\_ x Co x s )2As2- 

We study the vortex melting in Ba(Fei- x Co x )2As 2 (x=0.1) as a representative of 122- 
type iron-based superconductors [TTj . The parameter measured from a most recent 
experiment p] gives T c = 22K, \ ± = 160nm, £ ± = 2.44nm, H c c2 (T = 0) = 50T, 
H^(T = 0) = 70T, and s = 23A. The anisotropic parameter 7 falls from 2.0 to 1.5 
with the decrease of temperature. In our calculation, we observe that the melting lines 
determined in the framework of elastic theory only change slightly with the variation of 
7 (in the range of 1.5 ~ 2.0). So we set 7 = 2 in our calculations. 

We first calculate the melting line for the parallel fields with various Lindemann 
numbers. It is interesting to note in Fig. 1 that the vortex melting line for 
Ba(Fei_ x Co x )2As2 obtained in this paper is well consistent with the irreversible line 
measured experimentally [12] with Lindemann number c = 0.2. The irreversibility line 
in superconductor is usually regarded as the melting line. [30] 

Then we study the vortex melting for perpendicular fields. Due to the laking of 
the experimental data until now, three groups of Lindemann numbers are selected to 
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Figure 2. Vortex-lattice melting lines for Ba(Fei- x Co x )2As2 with magnetic field 
perpendicular to c-axis (without intrinsic layer pinning). From left to right, C x — 
0.05, C y = 0m,C v = 0.0375; C x = 0.1, Cj, = 0.1, C y = 0.075; C x = 0.2, C y = 
0.2, C y = 0.15. 
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Figure 3. Vortex-lattice melting lines for Ba(Fei_ x Co x )2As2 with magnetic field 
perpendicular to c-axis (with intrinsic layer pinning). 

calculate the melting lines (C x = 0.05, C y = 0.05, C y = 0.0375; C x = 0.1, C y = 0.1, C y = 
0.075; C x = 0.2, C y = 0.2, C y = 0.15). The obtained melting curves are shown in Fig. 2. 
The future experimental results can be inserted into this figure, and be compared with 
interpolation values obtained here. 

Setting C x = C y , we obtain two curves as shown in Figs. 2(a-c), as in Ref. [31]. As 
pointed out in Ref. [32], to interpret these two curves as two " melting lines " and thus 
reach to a conclusion of an intermediate phase is unphysical, since the elastic theory can 
give at best one melting line. To impose the same Lindemann number along the two 
directions is of no physical basis. In order to achieve a single melting line, one can tune 
the Lindemann numbers in the two directions. A good collapse of the melting lines in 
two directions can be achieved if setting the ratio c x /c y ~ 1.33. It is interesting to note 
that this ratio is very close to 1.37 observed in Ref. [33] using parameters of cuprate 
superconductors. 

The vortex melting for fields perpendicular to the c-axis is also influenced by the 
layer pinning. In order to study this intrinsic layer pinning effect, the matrix (15) 
is used to calculate the thermal fluctuations along two transverse directions. The 
melting lines for two groups of Lindemann number {C x = 0.1, C y = 0.1, C y = 0.075; 
C x = 0.2, C y = 0.2, C y = 0.15) are collected in Fig. 3. A good collapse can also be 
obtained by setting the same ration c x /c y ~ 1.33. The future experimental data can be 
also compared with these values. 
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Figure 4. Vortex-lattice melting lines for N d(Oi- x F x )FeAs with magnetic field 
parallel to c-axis. 



3.2. Vortex melting in Nd(0 1 _ x F x )FeAs. 

We now turn to another typical iron-based layered superconductor Nd(Oi- x F x )FeAs(llll- 
type) [ES El H51 HE|- We take parameters from PJ, T c = 49K, \ ± = 200nm, 
£i = 2.3nm, # C C 2 (T = 0) = 95T, F c a |(T = 0) = 220T, 7 2 = 25 . The separation 
between the FeAs planes is given by twice of the lattice parameter s = 2c ~ 17 A 
[DEI [22]. 

Similarly, we first calculate the melting lines for the fields parallel to 
shown in Fig. 4. Since the melting line has not been explicitly reported, the data at 
which the resistance reaches the small percentage (0.5%) of the normal-state resistance 
in Ref. |15J is regarded the approximated melting one here, which is also collected in 
Fig. 4. Interestingly, the theoretical melting line with c = 0.18 is consistent with the 
experimental data in a considerable wide temperature range. The discrepancy between 
experiment and theory becomes larger at high magnetic fields and low temperatures. 
A trend of temperature independent melting data is shown at high magnetic fields in 
experiments, implying that the dimensional crossover occurs at this regime. If the vortex 
system is of two dimensional nature at high magnetic fields, the melting line is almost 
independent of the temperature [30] . 

Next, We calculate the melting lines when the fields perpendicular to c-axis, which 
are shown in Fig. 5. For perpendicular fields, three different groups of Linderman 



numbers (C x = 0.1, C y = 0.1, C y = 0.075; C x = 0.2, C y 



0.2, C y 



0.15; C T = 0.3, a, 
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Figure 5. Vortex-lattice melting lines for N d{Oi- x F x )FeAs with magnetic field 
perpendicular to c-axis (without intrinsic layer pinning). From left to right, C x = 
0.1, C y = 0.1, Cy = 0.075; C x = 0.2, C y = 0.2, C y = 0.15; C x = 0.3, C y = 0.3, C y = 
0.225. 

0.3, C y = 0.225) are used to locate the melting lines. A good collapse of the melting 
lines in two directions can be also achieved if setting the ratio c x /c y ~ 1.33, nearly the 
same as that in 122 type. 

We also investigate the intrinsic layer pinning effect on the vortex melting in 
Nd(Oi- x F x )FeAs sample. The results are exhibited in Fig. 6. Surprisingly, we observe 
an intersection between two melting lines with two transverse directions for typical 
Lindemann numbers, which is not shown in the above 122 type. Obviously, a collapse 
of these two curves could not be achieved by setting any ration c x /c y . It is implied that 
the intermediate smectic phase [311 [321 [33] may exist in Nd{Oi_ x F x )FeAs. However, 
its identification is beyond the phenomenological Lindemann theory. It is interesting 
to note that, when the exponentially weak intrinsic pinning was taken into account for 
a weakly anisotropic Ba(Fe\- x Co x )2As2, the vortex melting behavior is only changed 
slightly, but for the strongly anisotropic Nd(Oi- x F x )FeAs, the intrinsic pinning would 
play a much more significant role. 

4. Summary 

In terms of the elastic theory, we have studied the thermal fluctuations in two typical 
iron-based superconductors Ba(F : ei_ x C 'o x ) 2 As 2 and Nd{Oi_ x F x )FeAs for magnetic 
fields parallel and perpendicular to the anisotropy axis. Using the parameters of 
these superconductors from experiments, the melting lines with various Lindemann 
numbers are composed. Interestingly, we can derive the melting line observed in 
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Figure 6. Vortex-lattice melting lines for N d(Oi- x F x )FeAs with magnetic field 
perpendicular to c-axis (with intrinsic layer pinning). 



both B a(F e\- x C o x )2As2 [12] and N d(0\- x F x )FeAs p3] for magnetic field parallel to 
the c— axis. Neglecting the layer pinning effect, it is found that thermal fluctuations 
normalized by vortex separations in the two transverse directions for parallel field are 
proportional to each other, similar to those observed in cuprate superconductors [321 133] . 
Interestingly, the intrinsic layer pinning play much more important role on the vortex 
melting in Nd{0\- x F x )FeAs than in Ba{Fe\- x Co x )2As2- More experimental works 
are motivated to confirm our prediction by comparing the melting dada with the 
interpolation one presented here. The present prediction may in turn provide a guide 
or reference to locate the melting line in future experiments. 
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